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ABSTRACT 


The  objective  of  this  thesis  was  to  develop  a  computer 
program  for  use  in  evaluating  data  acquired  from  the  Depart- 
ment of  Aeronautics  Cessna  310H  model  aircraft.   The  program 
was  designed  for  use  with  the  IBM  360  ^omputer,  in  conjunc- 
tion with  a  CALCOMP  plotter.  The  program  format  is  very 
similar  to  that  of  IBM's  Scientific  Subroutine  Package 
routines;  in  addition,  the  program  provides  the  student  with 
a  detailed  set  of  instructions  for  loading  the  data. 

Using  information  taken  from  the  engine  curves  for  the 
aircraft,  the  program  evaluates  input  data,  determines  the 
power  required  for  flight  at  altitude,  reduces  the  values 
to  equivalent  sea  level  values  and  outputs  a  graph  of  power 
versus  velocity  for  level  flight.  Lift  and  drag  coefficients 
are  calculated  as  functions  of  aircraft  weight,  velocity, 
and  power  required,  and  are  plotted  on  a  second  graph. 
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I.   INTRODUCTION 

The  laboratory  course  is  one  of  the  most  demanding  for 
the  graduate  student  because  of  the  extremely  time  consuming 
process  of  reducing  raw  data.   Often  the  student  may  find 
himself  neglecting  important  course  work  because  of  the 
large  amount  of  work  required  to  reduce  laboratory  data  by 
hand  to  meaningful  results.  This  is,  of  course,  an  unsat- 
isfactory situation,  but  is  one  which  seems  to  be  the  rule 
rather  than  the  exception  for  many  students. 

The  Department  of  Aeronautics  Cessna  310H  aircraft 
flying  laboratory,  which  supports  the  Flight  Evaluation 
Techniques  courses,  AE  3321  and  AE  3322,  at  the  Naval  Post- 
graduate school,  is  typical  of  such  time  consuming  labora- 
tories. With  the  incorporation  of  a  photo-panel  data 
acquisition  system  in  the  aircraft  came  the  capability  to 
collect  large  amounts  of  in-flight  data  which  require  large 
amounts  of  the  student's  time  for  evaluation. 

There  appeared  to  be  a  need  for  some  sort  of  procedure 
which  would  reduce  the  student  workload  in  this  laboratory 
by  relieving  him  of  the  burden  of  hand  reducing  the  large 
volume  of  raw  data  and  hand  plotting  the  results.   A  program 
which  would  require  only  that  the  student  assemble  the  data 
in  a  proper  sequence  and  input  it  to  a  computer,  which  in 
turn  would  evaluate  the  data  and  output  tabular  and  graphical 
displays  of  the  results,  would  be  extremely  beneficial  in 
lightening  the  student  workload.   Such  a  program  would  free 
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the  student,  enabling  him  to  spend  time  evaluating  his 
results  instead  of  tabulating  them.   With  this  need  in  mind 
a  computer  program  for  the  reduction  of  in-flight  data  for 
the  Cessna  310H  aircraft  was  written. 

CESSNA  was  designed  as  a  student  aid,  with  the  intent 
that  it  be  self-explanatory,  containing  all  information  and 
instructions  necessary  for  the  student  £o  run  the  program. 
Parts  II  and  III  of  this  thesis  may  be  used  as  a  user's 
manual  in  the  event  more  detailed  information  is  required, 
or  if  modification  and/or  additions  are  desired.  The  program 
is  written  in  FORTRAN  language  but  knowledge  of  the  language 
is  not  required,  other  than  the  ability  to  type  data  cards 
and  a  job  card  in  proper  format;  this  format  information  is 
readily  available  at  the  computer  center*,  The  card  deck  is 
Complete.   The  student  needs  only  to  add  his  data  cards  to 
the  end  of  the  deck  and  it  is  ready  to  run. 


II.   PROGRAM  ORGANIZATION  AND  OPERATION 

A.  GENERAL  ORGANIZATION 

As  shown  in  Fig.  1,  the  complete  program  consists  of 
two  primary  sections  and  two  subroutines: 

1.  Aircraft  Power  Requirements: 

2.  Lift  and  Drag: 

3 .  Interpolation  Subroutine: 
4«   Sorting  Subroutine. 

When  program  execution  begins,  aircraft  engine  data  are 
read  in  for  internal  use  in  the  program.  Following  this, 
student  data  concerning  aircraft  weight,  altitude,  velocity, 
RPM  and  manifold  pressure,  and  temperature,  are  read  in  for 
evaluation. 

CESSNA  uses  the  engine  data  and  the  student's  RPM  and 
manifold  pressure  data  to  determine  shaft  horsepower  required 
for  level  flight  at  the  operating  altitude  for  the  input 
velocity  range.   Horsepower,  along  with  velocity  and  weight, 
is  then  used  to  determine  lift  and  drag  coefficients  for 
the  aircraft.  Results  are  printed  out  in  tabular  form  and 
are  graphed  using  a  CALCOMP  plotter. 

B.  AIRCRAFT  POWER  REQUIREMENTS 

The  heart  of  the  program  is  the  set  of  engine  perfor- 
mance curves  shown  in  Fig.  2.   On  the  sea  level  side  of  the 
chart  the  curves  are  a  linearized  version  of  the  original 
curves  (Continental  Motors  Corporation,  195<3)  ,  and  are  a 
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good  approximation  to  the  originals.   In  the  recommended 
cruising  range  the  maximum  error  in  the  linearized  curves 
is  approximately  1  to  1.5$  in  the  mid-range  RPM's  and  mani- 
fold pressures  (or  a  maximum  of  approximately  2.5  brake 
horsepower),  and  much  lower  in  the  upper  and  lower  ranges, 
where  the  original  curves  are  very  nearly  linear.   The 
altitude  curves  are  a  copy  of  the  original  curves. 

The  curves  are  utilized  in  the  following  manner.  For 
each  RPM  line  in  the  sea  level  set  of  curves,  the  slope  of 
the  line,  maximum  manifold  pressure,  and  maximum  brake  horse- 
power have  been  determined.   For  the  altitude  set  of  curves 
each  RPM  line  has  had  slope  and  maximum  brake  horsepower 
determined.   A  question  may  arise  at  this  point  as  to  why, 
for  a  given  RPM,  the  maximum  brake  horsepower  is  not  the 
same  for  both  curves  and  why,  therefore,  one  horsepower 
value  is  not  sufficient  for  any  calculations  needed.   The 
answer  is  two-fold.  First,  in  linearizing  the  sea  level 
set  of  curves,  the  maximum  points  for  the  curves  changed, 
whereas  the  maximum  points  for  the  altitude  side  did  not; 
hence,  two  values.   Secondly,  since  points  on  the  RPM  lines 
are  determined  using  slopes  and  maximum  points,  it  is 
necessary  to  keep  both  values  so  that  calculations  for  both 
sets  of  curves  will  be  accurate. 

A  simplifying  assumption  was  made  concerning  the  mani- 
fold pressure  lines  on  the  altitude  side.  For  purposes  of 
this  program  the  lines  were  assumed  to  be  vertical.  This 
introduces  an  error  of  not  more  than  1.1$,  which  yields  at 
most  a- difference  of  two  brake  horsepower.   This  assumption 
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ENGINE  PERFORMANCE  CURVES 


i  ■? 


greatly  simplifies  the  method  by  which  manifold  pressure  is 
accounted  for.   Since  the  altitude  curves  are  graphed  as 
a  linear  function  of  the  density  ratio,  a  simple  conver- 
sion is  made  from  manifold  pressure  to  density  ratio  and 
this  converted  ratio  (called  MAPSIG)  is  used  for  calcula- 
tions, as  will  be  shown  later.  Table  I  shows  the  data  used 
for  the  MAPSIG  conversion.  This  is  a  purely  mathematical 
manipulation,  used  for  simplifying  the  equations  by  restating 
one  variable  in  terms  of  another.  This  converted  ratio  is 
related  to  the  critical  operating  altitude  and  is  there- 
fore differentiated  from  the  true  density  ratio  (SIGMA)  for 
the  actual  operating  altitude,  which  is  determined  in  the 
program  from  the  flight  altitude.  Table  II  shows  the  data 
used  for  determining  the  flight  altitude  density  ratio. 
Once  the  engine  curve  parameters  have  been  determined,  power 
requirements  are  calculated  for  a  single  engine. 

On  the  sea  level  side,  brake  horsepower  (SSLBHP)  is 
calculated  as  a  function  of  the  maximum  brake  horsepower 
(MAXBHP)  and  maximum  manifold  pressure  (MAXMAP)  for  the 
operating  RPM,  slope  of  the  RPM  line  (RSLOPE) ,  and  operating 
manifold  pressure  (MAP)  according  to  the  equation: 

SSLBHP  =  MAXBHP  -  RSLOPE  (MAXMAP  -  MAP)  (l) 

The  procedure  for  determining  a  value  for  brake  horsepower 
at  altitude  is  reasonably  straightforward.   The  critical 
altitude  (BHPCR)  for  a  given  manifold  pressure  is  that 
altitude  represented  by  the  intersection  of  the  operating 
RPM  and  manifold  pressure  lines.   The  converted  density 
ratio  (MAPSIG)  discussed  earlier  defines  this  altitude. 
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TABLE  I 

MAPSIG  CONVERSION  DATA 

MAP 

MAPSIG 

MAP 

MAPSIG 

MAP 

MAPSIG 

14. 

.545 

20. 

.730 

26. 

.905 

15. 

.520 

21. 

.760 

27. 

.935 

16. 

.615 

22. 

.790 

28. 

.965 

17. 

.640 

23. 

.820 

29. 

.990 

IB. 

.675 

24. 

.855 

19. 

.705 

25. 

.880 

ALT 


TABLE  II 
ALTITUDE/DENSITY  RATIO  CONVERSION  DATA 
SIGMA  ALT      SIGMA         ALT      SIGMA 


0. 

1.00000 

6000. 

.83586 

1000. 

.97106 

7000. 

.8IO64 

2000. 

.94277 

8000. 

.78601 

3000. 

.91512 

9000. 

.76196 

4000. 

.88808 

10000. 

.73848 

5000. 

.86167 

11000. 

.71555 

12000. 

.69317 

13000. 

.67133 

14000 c 

.65002 

15000. 

.62923 
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A  straight  line  is  drawn  connecting  this  point  with  the 
sea  level  brake  horsepower  value.   Then  a  vertical  line  is 
drawn  upward  from  the  operating  density  ratio  on  the  abscissa. 
If  this  line  intersects  the  SSLBHP-BHPCR  line,  then  that 
intersection  defines  the  brake  horsepower  required  for 
level  flight  at  the  operating  altitude.   If  the  vertical 
line  lies  to  the  right  of  the  critical  altitude  then  its 
intersection  with  the  operating  RPM  line  defines  the 
required  brake  horsepower.  Figure  3  shows  an  example  of 
engine  curve  utilization. 

Critical  brake  horsepower  (BHPCR)  is  calculated  as  a 
function  of  maximum  brake  horsepower  for  the  operating  RPM 
(MABHP),  the  slope  of  the  operating  RPM  line  (ARSLOP)  and 
the  converted  density  ratio  (MAPSIG)  according  to  the 
equation: 

BHPCR  =  MABHP  -  ARSLOP   (1.0  -  MAPSIG)       (2) 

The  slope  of  the  line  connecting  sea  level  and  critical 
brake  horsepowers  (LSLOPE)  is  also  a  function  of  the  con- 
verted density  ratio: 

LSLOPE  =  (BHPCR  -  SSLBHP)  /  (1.0  -  MAPSIG)      (3) 

The  operating  altitude  is  compared  with  the  critical  altitude. 
If  the  operating  altitude  is  the  lower  one,  brake  horsepower 
is  calculated  on  the  operating  RPM  line  as: 

BHPALT  =  SSLBHP  +  LSLOPE  (1.0  -  SIGMA)  (4) 
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If  the  operating  altitude  is  above  the  critical  altitude 
then  brake  horsepower  is  calculated  on  the  operating  RPM 
line  as: 

BHPALT  =  MABHP  -  ARSLOP  (1.0  -  SIGMA)        (5) 

It  should  be  noted  at  this  point  that  the  actual  flight 
data  are  based  on  thrust  horsepower  rather  than  brake  horse- 
power.  In  the  classroom  the  student  takes  an  ideal  brake 
horsepower,  applies  the  necessary  corrections  and  then  mul- 
tiplies that  value  by  propeller  efficiency  to  get  the  actual 
thrust  horsepower  available.  This  program  does  exactly  the 
opposite.   From  in-flight  RPM  and  manifold  pressure  it  cal- 
culates the  actual  thrust  horsepower  required.   This  is 
corrected,  and  then  converted  to  shaft  horsepower  (it  is 
assumed  that  shaft  and  brake  horsepower  are  equal)  by  divid- 
ing by  propeller  efficiency. 

Since  the  engine  performance  curves  are  plotted  in 
terms  of  brake  horsepower,  the  nomenclature  of  the  program 
up  to  this  point  has  been  written  in  terms  of  brake  horse- 
power. However,  since  the  horsepower  value  just  found  using 
equation  (5)  is  actually  a  thrust  horsepower,  from  here  on 
it  will  be  referred  to  as  a  thrust  horsepower,  until  it  is 
converted  to  shaft  horsepower. 

Once  the  basic  thrust  horsepower  value  has  been  deter- 
mined, several  correction  factors  are  applied,  the  first 
of  which  is  for  temperature  variation.   The  engine  curves 
are  plotted  for  standard  temperature;  a  correction  factor 
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Figure  3 
ENGINE  CURVE  UTILIZATION 
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is  supplied  with  them.   First  the  temperature  is  changed 
from  degrees  Fahrenheit  (TEMPF)  to  degrees  Centigrade 
(DEGC)  and  the  difference  between  operating  and  standard 
temperatures  (TDIF)  is  determined: 

DEGC  =  5.0  (TEMPF  -  32.0)  /  9.0  (6) 

TDIF  =  15.0  -  DEGC  \  (7) 

Then,  using  the  correction  factor,  the  single  engine  thrust 
horsepower  corrected  for  temperature  variation  (CBHP)  is: 

CBHP  =  BHPALT  (1.0  -0.0016667  (TDIF))       (g) 

The  total  thrust  horsepower  (TBHPA)  is  found  by- 
multiplying  the  temperature  corrected  value  by  the  number 
of  engines: 

TBHPA  =  CBHP  (NOENG)  (9) 

Total  shaft  horsepower  required  at  altitude  is  found 
by  dividing  thrust  horsepower  by  propeller  efficiency: 

TSHPA  =  TBHPA  /  PE  (10) 

Propeller  efficiency  has  been  arbitrarily  chosen  as  0.75 
and  is  entered  in  the  program  using  a  DATA  statement  for 
ease  of  access,  should  subsequent  evaluation  dictate  a 
change  in  the  value. 

The  total  shaft  horsepower  is  good  for  one  altitude 
and  one  temperature.   To  make  this  information  more  useful 
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it    is   converted  to   an   equivalent    shaft   horsepower   (SSLPWR) 
which  may  be  used  to  find  horsepower  at   any  altitude  for 
the  given  temperature: 


SSLPWR   =  TSHPA        / SIGMA  (11) 

The   relationship  between  velocity  and  power  and  their   equiv- 
alent   counterparts   is  most    easily   shown  by  the  following 
sequence   of   equations: 

*ssl  =  *>o      V" 

and 

lalt  =  *    1°      V 
also 

or   =  /y>0  ;  /=>  =  cr  f>Q 

so 

qalt=i^0crV2 
Define   equivalent   velocity:      V      =  fij~     V 

then 

*alt   =  *       /°o        Ve2 

P   =  TV   =   DV 

=  i  S  cd  ^o  V2  V  =  i  S  cd  p   V3 

=  *  S  Cd  ^  cr   V3 
Define   equivalent   power:      P      =  P  ^CT 

p  /F"        =  i  s  cd/o  Jk  v3    =  i  s  cd  ^  (  y^r  v)3 
=  i  s  c       ^       v  3 

d       'o  e 
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Once  equivalent  power  has  been  calculated,  the  remaining 
step  in  the  power  requirements  section  is  to  change  velocity 
to  equivalent  velocity  (EQVEL): 

EQVEL  =  VEL  J  SIGMA  (12) 

C.    LIFT  AND  DRAG 

Lift  is  a  function  of  density,  velocity,  characteristic 
area  and  the  lift  coefficient.   In  order  to  avoid  the 
necessity  of  inputing  or  calculating  a  density  for  the 
flight  altitude,  the  equations  for  determining  lift  and 
drag  coefficients  are  based  on  equivalent  velocity,  reducing 
density  to  a  constant  (sea  level  value). 
Lift  =  weight  =  ^  /o  V2  S  C1 

=  i  /°0  Or    V2  S  C1 

=  i  f>0  (^v)2  s  c1 

=  *  /%   Ve2  S  Cl 

For  the  calculations,  equivalent  velocity  is  changed  from 
knots  to  feet  per  second  (VEFPS): 

VEFPS  =  EQVEL  (6000.)  /  3600.  (13) 

Then  the  lift  coefficient  (CLIFT)  is  found  as  a  function 
of  weight  (WT)  and  equivalent  velocity: 

CLIFT  =  2  (WT)  /  (.0023769)  (S)  (VEFPS)2     (14) 
The  characteristic  area,  S,  is  defined  by  a  data  state- 
ment at  the  beginning  of  the  program  for  ease  of  access. 
There  are  two  values  given  for  S  (NASA  TN-D-623S,  1971). 
In  this  program  200  square  feet  is  used,  which  includes  tip 
tanks.   Without  them  the  wing  area  is  179  square  feet. 
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Drag  is  a  function  of  density,  velocity,  characteristic 
area  and  the  drag  coefficient.   It  is  also  a  function  of 
velocity  and  power: 

Drag  =  P  /  V  =  Pe  /  Ve 

=  i  p   V2  S  Cd 
or  equivalently,    P     =  -§■  p    S  C,  V     V 

=  */^oSCdVe3 

For  the  calculations,  equivalent  power  is  changed  from 
shaft  horsepower  to  foot-pounds  per  second  (PEFPPS): 

PEFPPS  =  SSLPWR  (550.)  (15) 

Then  the  drag  coefficient  (CDRAG)  is  found  as  a  function 
of  equivalent  power  and  equivalent  velocity: 

CDRAG  =  2  (PEFPPS)  /  (.0023769)  (VEFPS)3    (l6) 
Once  all  computations  have  been  completed,  a  table  of 
velocity,  manifold  pressure,  brake  horsepower  at  altitude, 
equivalent  velocity,  equivalent  shaft  horsepower,  lift 
coefficient  and  drag  coefficient  is  printed  out.   Then  two 
calls  are  made  to  DRAW,  one  of  IBM's  plotting  subroutines 
used  in  conjunction  with  a  CALCOMP  plotter.   The  first  call 
is  for  a  graph  of  equivalent  power  versus  equivalent  velo- 
city; the  second  is  for  a  graph  of  lift  coefficient  versus 
drag  coefficient. 

D.    INTERPOLATION  SUBROUTINE 

As  stated  earlier,  the  heart  of  the  program  is  the 
data  obtained  from  engine  curves  for  the  Cessna.   Originally 
it  was  decided  that  these  data  would  be  used  to  formulate 
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equations  with  RPM  as  the  independent  variable.  Then, 
given  an  operating  RPM,  all  necessary  information  could  be 
obtained  using  the  equations,  eliminating  the  need  for 
table  searches  and  interpolation.   Several  subroutines 
from  IBM's  Scientific  Subroutine  Package,  including  LSQPL2, 
CHBFT,  and  SPLIN,  were  utilized  in  an  effort  to  formulate 
equations  which  would  accurately  represent  the  engine  data. 
However,  due  to  problems  with  convergence,  extended  pre- 
cision requirements  for  curve  coefficients,  and  decreasing 
accuracy  with  increasing  curve  order,  all  of  these  routines 
proved  unsatisfactory  in  producing  the  required  equations. 
It  was  therefore  decided  to  write  a  simple  subroutine  which 
would  interpolate  in  a  table  of  data  from  the  engine  curves. 

Subroutine  LINT  performs  linear  interpolation  or 
extrapolation.   Inputs  to  the  subroutine  are  two  columns 
of  data  and  an  independent  variable  for  which  LINT  will 
interpolate  the  dependent  variable.  The  data  column 
corresponding  to  the  independent  variable  must  be  in 
monotonically  increasing  order  but  may  be  arbitrarily 
spaced.   Since  the  engine  curves  are  for  a  single  engine, 
flight  with  two  engines  allows  operation  of  each  engine  at 
values  below  those  shown  on  the  engine  performance  curves, 
but  not  above.   Subroutine  LINT  was  therefore  written  to 
extrapolate  below  the  limits  of  the  curves  but  not  above. 
In  the  event  that  the  independent  variable  falls  above 
the  range  of  data  given,  an  error  message  so  stating  this 
is  printed  out. 
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The  method  used  by  LINT  is  very  simple.   First  a  check 
is  made  to  see  if  the  variable  is  below  the  range  for  inter- 
polation; if  this  is  so,  an  extrapolation  is  performed. 
Otherwise  the  column  is  searched  until  bracketing  values 
are  found,  at  which  time  an  interpolation  is  performed.   If 
a  bracketing  pair  cannot  be  found,  the  variable  is  above 
the  data  range  and  the  error  message  is  printed. 

E.    SORTING  SUBROUTINE 

To  make  the  program  easier  to  use,  it  was  written  so 
that  the  student's  data  may  be  loaded  in  any  order.   It  is 
not  necessary,  for  example,  to  put  the  data  in  order  of 
increasing  velocity.   Subroutine  SORT  does  this  ordering 
during  program  execution  so  that  the  printed  output  data 
is  ordered  in  terms  of  increasing  velocity.   SORT  is 
written  to  evaluate  one  column  of  data,  and  to  put  two 
associated  columns  of  data  in  the  same  chronological  order 
as  the  one  being  evaluated. 

The  method  used  by  SORT  is  simple.   A  logical  flag 
is  set  equal  to  FALSE.   Then  sequential  pairs  of  values 
in  the  first  data  column  are  compared.   If  the  second 
value  in  any  one  of  these  pairs  is  less  than  the  first 
value  in  the  pair,  the  positions  of  the  two  values,  and 
of  the  associated  values  in  the  other  two  columns,  are 
switched  and  the  logical  flag  is  set  equal  to  TRUE. 
When  the  entire  column  has  been  searched,  the  flag  is 
checked;  if  it  has  been  reset  (to  TRUE),  this  indicates 
the  possibility  that  there  may  be  more  switches  necessary 
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so  the  flag  is  reset  to  FALSE  and  the  procedure  is  repeated. 
When  a  complete  column  search  is  made  without  resetting 
the  flag,  the  data  are  in  proper  order  and  control  returns 
to  the  main  program.  Like  the  main  program,  the  subroutine 
is  limited  to  100  data  points  per  column. 

III.   PROGRAM  USAGS, 

A.    DATA  LOADING 

The  first  data  to  be  loaded  following  execution  of 
the  program  are  the  engine  data,  consisting  of  the  slope, 
maximum  manifold  pressure  and  maximum  brake  horsepower 
(again,  brake  horsepower  values  are  used  here  because  the 
engine  data  comes  from  curves  which  are  plotted  in  terms 
of  brake  horsepower)  for  sea  level  performance,  and  slope 
and  maximum  brake  horsepower  for  altitude  performance,  for 
each  of  nine  sets  of  RPM  lines.  These  data  are  shown  in 
Table  III.   Following  the  permanent  data  are  the  student's 
in-flight  data. 

The  order  in  which  in-flight  data  is  loaded  varies. 
The  program  gives  the  option  of  constant  or  variable  RPM, 
and  the  input  order  depends  on  which  option  is  selected. 
To  simplify  the  student  data  loading,  input  data  formats 
for  both  options  are  included  in  the  program.   These  formats 
are  similar  to  those  used  by  James  Melsa  in  his  book  of 
computer  programs  for  linear  control  theory.   They  give 
card  by  card  instructions  for  loading  the  data,  including 
card  number,  column  numbers,  data  description,  and  input 
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TABLE   III 

ENGINE   PERFORMANCE   DATA 


RPM 

RSLOPE 
(BHP/INCH) 

MAXMAP 
(INCHES   HG) 

MAXBHP 
(HP) 

ARSLOP 
(BHP/INCH) 

MABHP 
(HP) 

1300 

6.45 

29.75 

150.0 

173.91 

150.0 

1900 

7.16 

29.70 

164.5 

136.67 

163.0 

2000 

7.37 

29.65 

175.0 

200.00 

175.0 

2100 

3.17 

29.60 

193.0 

216.67 

191.0 

2200 

3.72 

29.55 

203.5 

235.16 

207.5 

2300 

9.24 

29.50 

223.5 

254.94 

223.5 

2400 

10.00 

29.40 

239.0 

267.03 

237.5 

2500 

10.56 

29.40 

243.0 

231.32 

249.0 

2625 

11.09 

29.35 

260.0 

294.50 

260.0 

TABLE   IV 
EXAMPLE   PROBLEM  DATA 


RPM  MAP 

(INCHES   HG) 


2100 
2100 
2100 
2100 
2200 
2300 
2300 
2200 
2100 
2100 
2100 
2100 
2450 


17.5 
19.0 

19.75 

20.3 

22.0 

23e6 

23.6 
22.0 
20.1 
19.2 
13.6 

17.5 
23.0 


VEL 
(KNOTS) 


RPM  MAP  VEL 

(INCHES   HG)         (KNOTS) 


35.0 
100.0 
115.0 
130.0 
145.0 
160.0 
160.0 
145.0 
130.0 
115.0 
100.0 

35.0 
164.0 


2450 

2300 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
2300 
2300 
2300 


20.0 

13.0 
13.0 
13.0 
15.0 
16.0 
18.0 
13.5 
19.0 
21.0 
21.5 
23.0 
24,0 


151.5 
135.0 

119.5 
103.0 

37.5 
30.0 
39.0 
111.5 
122.0 
130.5 
144.5 
152.5 
153.5 
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format.  These  formats  are  presented  in  the  program  in 
Appendix  A. 

Following  the  input  of  constant  flight  parameters, 
such  as  altitude  and  weight,  velocities  and  manifold  pres- 
sures are  read  in.   If  RPM  is  variable,  its  values  are 
read  in;  otherwise  the  constant  RPM  value  is  read  in  and 
duplicated  for  each  set  of  velocity  and  pressure.  This 
duplication  is  done  for  use  in  the  sorting  subroutine, 
which  requires  a  column  of  RPM  data,  and  for  the  printed 
output,  which  includes  a  column  for  RPM.   Once  all  the 
data  are  read  in  and  sorted,  the  engine  curve  performance 
parameters  are  interpolated  for  each  RPM,  if  RPM  is  vari- 
able; if  it  is  constant  the  parameters  are  evaluated  only 
once,  saving  unnecessary  computer  time.  Following  this 
evaluation,  the  power,  lift  and  drag  requirements  are 
calculated. 

B.    GRAPHICAL  OUTPUT 

The  power  requirements  and  lift  and  drag  coefficients 
for  the  aircraft  are  graphed  on  the  CALCOMP  plotter,  using 
IBM's  DRAW  subroutine,,   Each  graph  is  eight  inches  by 
ten  inches,  with  a  title  below*   The  abscissa  and  ordinate 
are  automatically  scaled  by  the  subroutine  and  the  scale 
information  is  printed  below  the  graph.   Data  points  are 
plotted  on  the  graph  as  triangles.   This  was  done  for 
several  reasons.   If  the  subroutine  is  instructed  to  plot 
a  continuous  curve,  it  merely  connects  each  point  to  the 
next  with  a  straight  line.   This  will  produce  a  very 
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erratic  curve  in  the  case  of  data  with  any  degree  of  spread, 
and  where  there  are  only  a  few  data  points  to  plot.   The 
use  of  triangles  allows  the  student  to  draw  a  smooth  curve 
through  the  data.   It  also  allows  the  data  to  be  loaded  in 
any  order,  rather  than  forcing  the  student,  for  example,  to 
arrange  his  data  in  increasing  order.   This  last  point  is 
admittedly  academic,  since  addition  ofVshe  sorting  sub- 
routine eliminates  any  plotting  problems  which  might  have 
occurred  from  unordered  data. 

IV.   EXAMPLE  PROBLEM 

Two  actual  sets  of  student  data  were  combined  for  the 
example  problem.  RPM  was  variable.  Table  IV  shows  the  RPM 
values  and  associated  manifold  pressures  and  velocities. 
Temperature  varied  slightly  during  the  test  flights;  an 
average  value  was  calculated  for  use  in  the  program.   Air- 
craft and  flight  parameters  are  shown  in  Table  V.   The 
order  of  the  data  deck  is  shown  in  Appendix  B. 

Results  of  the  program  compared  very  favorably  with 
hand  calculated  data.   Table  VI  shows  the  tabulated  output 
from  the  program.  There  was  some  spread  in  the  data,  which 
is  easily  visible  in  the  sample  graphical  output  shown  in 
Appendix  C  «, 
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TABLE  V 
EXAMPLE  PROBLEM  -  AIRCRAFT  FLIGHT  PARAMETERS 

Number  of  engines  operating:   2 
Flight  altitude  (feet):   5000 
Temperature  (degrees  F):   76.55 
Aircraft  weight  (pounds):   4500 
Number  of  data  points:   26 
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V.   PROJECT  EVALUATION 

It  is  felt  that  the  program  has  met  the  objectives 
previously  set  forth.   A  simple,  easy-to-use  program  has 
been  written  which  evaluates  aircraft  in-flight  data  and 
produces  graphs  of  power  versus  velocity,  and  lift  coeffi- 
cient versus  drag  coefficient.   It  yields  results  which  are 
comparable  in  accuracy  to  results  calculated  by  hand. 
Additionally,  there  is  a  considerable  savings  in  man-hours 
required  by  the  student  to  evaluate  his  data,  which  was  the 
original  impetus  for  this  thesis. 

VI.   AREAS  FOR  PROGRAM  EXPANSION 

There  are  several  possibilities  for  future  expansion 
of  the  program.  Temperature  variation,  though  small,  may 
have  a  noticeable  effect  on  performance  figures,  and  a  pro- 
gram option  allowing  variable  temperatures  would  be  bene- 
ficial.  The  same  is  true  for  an  option  which  would  allow 
altitude  to  vary.  Test  flights  are  subject  to  weather,  and 
clouds,  for  example,  do  not  always  allow  a  flight  to  be 
performed  at  a  single  altitude.   Another  variable  factor  is 
weight.  Fuel  accounts  for  nearly  16  percent  of  the  aircraft 
maximum  gross  weight,  and  on  a  long  flight  this  could 
appreciably  alter  performance  data.   An  option  which  uti- 
lizes fuel  readings  might  well  be  a  worthwhile  addition. 
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Another  section  which  might  be  added  to  the  program  is 
one  which  would  plot  rate  of  climb  curves  for  the  aircraft. 
This  could  be  a  simple  subroutine  which  would  require  inputs 
of  time  to  climb  through  a  given  altitude  change  over  a 
range  of  velocities.   An  easy  conversion  to  feet  per  minute 
could  be  made,  and  the  results  plotted. 

Other  future  possibilities  include  the  area  of  dynamic 
stability  and  control  information,  such  as  roll  rates,  which 
are  covered  in  the  second  half  of  the  flight  evaluation 
techniques  course  sequence. 
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APPENDIX    A 
COMPUTER  PROGRAM 
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APPENDIX  B 
SAMPLE  INPUT  DATA  DECK 


ir\    u"n 


• 

• 

• 

O 

LfN 

O 

O 

O 

ir\ 

ON 

-4- 

• 

• 

• 

O 

O 

o 

-4- 

H 

-4 

CM 

to 

H 

CM 

H 

Cf\ 

H 

H 

H 

CM 

H 

CM 

CM 

CM 

CM 

IT\  •         •         • 

•       •       •  vO  O     O     O 

o    u^   o  •••  o    o    o 

\C>     cn    oa  r^tOH  mc^iH 

HHH  CM      H     CM  CMCMCM 


IT\  •         •         • 

•       •       •  \D  OOO 

O     H     CM  •••  Ou^O 

sOlt\CM  o^OO  oa-4H 

i-H     H      H  CM      CM      H  CMCMCM 


ir\  •       •       o 

•       ••  O            u-\  OOO 

ua^-H  •••  Ou^O 

-4    v£>      H  CMr^tO  CM-4-H 

HvHH  CMCMH  CMCMCM 


•  to     u~\  OOO 

O  •  •  •••  OOO 

n  iA  O  OC^-tO  H     iH     H 

H  tO  tO  CMrHH  CMCMCM 


UTN 

• 

e 

• 

• 

• 

C^ 

vO 

O 

o 

O 

ur\ 

O 

• 

• 

o 

o 

O 

o 

o 

H 

o 

O 

o 

•to 

MD 

H 

H 

H 

H 

H 

-co 

H 

H 

H 

CM 

CM 

CM 

ITS  UN  •  •  •  • 

•  •  e       LT\  •  CM  OOOO 

\0  Ol^«*0'»»»«OOOO 

O  OHO-LTNONON^^-cfiHHHr^ 

tHH-COHiHHHCMCMCMCMCM 


©  •  Lf\  o  e  •  • 

OO  •       •       •     ia    H  OOOO 

OO  ••O'C0CMe»c.oOOO 

•     Oua  vOuAr^vOu^l^O'COoAr-IHrHr^v 

CMUA-4-OCMtOHHHiHCMTHCMCMCMCMCM 


43 


APPENDIX   C 
SJAMPLE  PROGRAM  GRAPHICAL  OUTPtJT 
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APPENDIX  D 
PROGRAM  GLOSSARY 

ALT  -  (input)  Altitude  (feet) 

ARSLOP  -  Slope  of  any  altitude  performance  RPM  line 

BHPALT  -  Single  engine  uncorrected  thrust  horsepower  at 
altitude 

BHPCR  -  Critical  horsepower  at  altitude  for  any  given  RPM 
and  MAP 

CI  -  Conversion  factor  for  changing  velocity  to  feet  per 
second 

C2  -  Constant  for  use  in  determining  drag  coefficient 

C3  -  Constant  for  use  in  determining  lift  coefficient 

CBHP  -  BHPALT  corrected  for  temperature  variation 

CDRAG  -  Drag  coefficient 

CLIFT  -  Lift  coefficient 

DEGC  -  Temperature,  degrees  Centigrade 

DUMMY  -  Dummy  variable  used  for  temporary  storage 

EQVEL  -  Sea  level  equivalent  velocity  (knots) 

IALT  -  Integer  version  of  ALT 

IPTS  -  Integer  version  of  POINTS 

ITB  -  Vector  of  constants  for  use  in  DRAW  subroutine 

IWT  -  Integer  version  of  wT 

LINT  -  Linear  interpolation  subroutine 

LSLOPE  -  Slope  of  the  line  connecting  sea  level  and  altitude 
values  of  BHP  (SSLBHP  and  BHPCR)  for  any  given  RPM 
and  MAP 

MABHP  -  Maximum  horsepower  for  any  given  altitude  performance 
RPM  line 

MAP  -  (input)  Manifold  pressure  (inches  Hg) 
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MAPSIG  -  Performance  curve  equivalent  value  of  density  ratio 
for  any  given  value  of  MAP 

MAXBHP  -  Maximum  horsepower  for  any  given  sea  level  perfor- 
mance RPM  line 

MAXMAP  -  Maximum  manifold  pressure  for  any  given  sea  level 
performance  RPM  line 

NOENG  -  (input)  Number  of  engines  operating 

ONERPM  -  (input)  Const ant /variable  RPM  flag 

PE  -  Propeller  efficiency 

PEFPPS  -  Sea  level  equivalent  power  (foot  pounds  per  second) 

POINTS  -  (input)  Number  of  data  points 

RA  -  Vector  of  altitudes  used  for  converting  altitude  to 
density  ratio 

RC  -  Matrix  of  engine  performance  curve  data 

RC1-RC6  -  Vector  equivalents  of  data  contained  in  RC 

RMAP  -  Vector  of  manifold  pressures  used  for  MAPSIG  con- 
version 

RPM  -  (input)  Engine  revolutions  per  minute 

RS  -  Vector  of  density  ratios  used  for  converting  altitude 
to  density  ratio 

RSIG  -  Vector  of  density  ratios  used  for  MAPSIG  conversion 

RSLOPE  -  Slope  of  any  sea  level  performance  RPM  line 

RTB1  -  Vector  of  data  used  in  DRAW  subroutine 

RTB2  -  Vector  of  data  used  in  DRAW  subroutine 

S  -  Aircraft  characteristic  area  (square  feet) 

SIGMA  -  Density  ratio  for  a  given  altitude 

SORT  -  Sorting  subroutine 

SSLBHP  -  Sea  level  power  for  any  given  RPM  and  MAP 

SSLPW  -  Sea  level  equivalent  total  shaft  horsepower 

TBHPA  -  Total  horsepower  at  altitude  corrected  for  tempera- 
ture and  number  of  engines 
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TDIF  -  Difference  between  standard  temperature  and  tem- 
perature at  flight  altitude 

TEMPF  -  (input)  Temperature  at  altitude  (degrees  Fahrenheit) 

TITLE1  -  Graph  title  information  for  DRAWP  subroutine 

TITLE2  -  Graph  title  information  for  DRAW  subroutine 

TSHPA  -  TBHPA  corrected  for  propeller  efficiency 

VEFPS  -  Sea  level  equivalent  velocity  (feet  per  second) 

VEL  -  (input)  Velocity  (knots) 

"WT  -  (input)  Aircraft  weight  (pounds) 
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